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7.8 Transient motion in a two-layered sea 
7.8.1 Shallow water approximation 




V5' 



y' =r-crr- 



Figure 7.8.1: A two-layered sea of contant depth. 

Let h be the mean depth and U, V be the depth-integrated fluxes in the upper layer, and 

h', U', V be the corresponding quantities in the lower layer. Let ( and (' be the vertical 
displacements of the free surface and interface respectively. Then 

(U,V)= f (u,v)dz, (U',V')= f k (u',v')dz, (7.8.1) 

J — h J — h—h' 

For shallow layers, we have hydrostatic pressure in the upper layer: 

P = Pg(( - z) (7.8.2) 

On the interface, z = —h + (', 

V = pg(C + h-C) (7.8.3) 
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In order that pressure is continuous on the interface we must have in the lower layer 

V = P9(( + h-C) + p'g(-h + C'-z) (7-8.4) 

Thus ^ g 

dx ^ dx (7 8 5) 

dx ^ \dx dx J ^ 9 dx ^ dx ^ ^ 9 dx 
The depth-integrated momentum equations in the upper layer are 

dU 1 dv dC t s 

7j7-f V = —T = -9hjr + - (7.8.7) 
dt pdx dx p 

dt pdy dy p 

9(C-C) 9U dV n , nnn , 
~dt + ~d^ + ~dy~ = (7A9) 

after using (7.8.2). Ignoring shear stresses on the interface and at the bottom, we have, for 

the lower layer, 



dU' 
dt 


-fv' = - 


1 dp' 
p' dx 


-gh -TT- 

p dx 


- tgh 

dx 


(7.8.10) 


dV 
dt 


+ fU' = - 


1 dp' 
p' dy 


-gh —tt- 

p 1 dy 


- egh — 

dy 


(7.8.11) 




d(' 


dU' 


f =o 

dy 




(7.8.12) 




dt 


dx 





where (7.8.2) has been used. The density contrast 

t = P —^- (7.8.13) 
P 

is usually small in oceans. 

There six equations for six unknowns : (, U, V, (', U', and V. As in the last section we can 
eliminate some unknowns and reduce the number of equations. For example, by eliminating 
U, V in the upper layer, we can get 

d ( d \ A u n ah^ dc - 1 9 ( dT * + dr y) f ( dT y dr ") (7 8 1^ 

Similarly by eliminating U', V in the lower layer we get, 

These two govern the coupled behavior of ( and (', in terms of which the no-flux boundary 
conditons on the velocity can also be expressed. 

We shall however introduce an alternate approach which has some advantages. 
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7.8.2 Transformation into normal form 

To avoid solving six simultaneous equations at once, we seek a linear combination so that 
the governing equations look like those for a one-layer fluid with just three unknowns. This 
is called the normal form. Take a(7.8.7) + (7.8.10), we get, 



Similarly a(7.8.8) + (7.8.11) gives 

!<„„ + n + /( at / + £0 = - 9 /»| (< + 0C + £ <) + f (7.8.17) 



Finally a(7.8.9) + (7.8.12) gives, 

-« - a£ + C') + ^(aU + U') + ^(aV + V') = Q (7.8.18) 

from the continuity equations. 
Now if we take 

and 

aC - aC + C = a( + (1 - a)C' = C (7.8.20) 
Then (7.8.16)-(7.8.18) become 

dU - dl f s 

-- f v = - 9 p h S + a ^ (7 . 8 . 21) 

8V - dl f s 

^r + fU = -gPh^- + a^- (7.8.22) 
at ay p 

dC dU dV 

and take on the appearance of the governing equations for a single layer fluid. 
For this reduction to be possible, conditions (7.8.19) and (7.8.20) imply that 

< + + j< = P K - oC + CO (7.8.24) 

Equating the coefficients of ( and £' separately, we get 

p a = a+ 1 lP- (7.8.25) 
Zip' 



4 



P(-a + l) 



h'e 
~~h 



(7.8.26) 



A quadratic equation is found for [3 by eliminating a 

h'e h'p\ h'e 



^-/M 1 + - + ^) + t = 



h hp' 



i.e., 



'h' 



P- 7T + 1 l3 + 



h!e 
~~h 



There are two solutions for (3 



Pi 

fa 



h + h' 
2h 



l±Wl- 



Aehh' 
(h + h') 



(7.8.27) 



(7.8.28) 



to be refered to as two normal modes. The solutions corresonding to Pi and P2 will be refered 
to as Mode 1 and Mode 2, respectively. 

Let us see what the corresponding a's are. Using (7.8.26), we have 



h'e 

a=l ~w 



(7.8.29) 



Since from (7.8.27), 



it follows that 



P1P2 



a = 1 



hfe 
~h 

P1P2 
P ' 



Therefore, 

ai = 1 - P 2 , a 2 — 1 - Pi 
Let us summarize the normal form equations for Mode k with k — 1, 2. Defining, 

U k = a k U + U', V k = a k V + V 
Ck = a k ( + (l-a k )(' 
(f x )k = a fc rf , (f y ) k = a k Ty. 



then 



at ox p 

at ay p 

da + dU k + dV k = Q 
dt dx dy 



(7.8.30) 

(7.8.31) 
(7.8.32) 
(7.8.33) 

(7.8.34) 
(7.8.35) 
(7.8.36) 
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Note that for mode k the effective gravity is g(3k, and the effective forcing stresses are 
and dkTy . Solution for the set of three equations is obviouly a simpler task. For a general 
wind stress field (r,f , Ty), we first solve for (Ui, Vi, (j) and (U2V2, C2) separately from (7.8.34) 
to (7.8.36). Afterwards (U,V,Q and ([/', V , (,') can be solved from the algebraic equations 
(7.8.31) and (7.8.32). For example, From 



Ci = aiC + (1 - ai)C', C2 = a 2 C + (1 - a 2 )C 
we get for the sea surface, 



c = 



Ci 


1 


— Oi 




C2 


1 


- «2 


(1 - a 2 )Ci - (1 - ai)C2 


dl 


1 


- Oi 


ai — a 2 


a 2 


1 


— a 2 





and for the interface 



c 





ai Ci 








a 2 C2 




-a 2 Ci + aiC 2 


ax 1 — ai 


ai — a 2 


a 2 1 — a 2 





Similarly, from 

we get, for the upper layer 



U 1 = a 1 U + U', U 2 = a 2 U + U' 



U = 



Ui 


1 








1 


Ui 


-u 2 




1 


ai 


- a 2 


a 2 


1 







and for the lower layer 



U' = 



ai Ui 
a 2 U 2 



CLi 1 

a 2 1 



-a 2 Ui + aiU 2 
ai — a 2 



The formulas for V and V are similar, with U's replaced by Vs. 

7.8.3 Normal form for small density difference 

From here we shall only be interested in small density differences, then 

h' eh' 



(3 2 = 



h h + h 
eh' 



Ot?) 



h + h 



7 + 0(6 2 ) 



(7.8.37) 



(7.8.38) 



(7.8.39) 



(7.8.40) 



(7.8.41) 



(7.8.42) 



(7.8.43) 
(7.8.44) 
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or 

eh' 

01 = 1 " aTa 7 + ° (e2) ' (7 ' 8 ' 45) 

Taking only the leading order terms, we have 

a ± = 1, a 2 = ~. (7.8.47) 

From (7.8.43) 



We also get from (7.8.41) and (7.8.42) 



c- (i+Av ^r i)& ^. 

v- v ,_ {h'/h)V 1 + V 2 

v ~l+h'/h' v - l + h'/h (7 - 8 - 52j 

For Mode 1 we get from (7.8.31), (7.8.32) and (7.8.33), 

jj x = U + U' Vi — V + V' 

Ci = C, (r,)i = rf (tv)i = (7-8-53) 
Let H = h + h' denote the total depth, the normal mode equations are 

§ +/M f + f ( ,, 55) 

<9t dx dy 

In view of (7.8.53), the two layers move together as if a single layer of homogeneous fluid. 
This is called the Barotropic (surface) mode. 
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For Mode 2 we recall from (7.8.34), 



a2 = 4 + 0(£ ), = ^ (7.8.57) 
Then from (7.8.31) to (7.8.33), 

U 2 = -^-U + U', v 2 = -^v + v 

h h 

and 

& = ^ (7.8.59) 
from (7.8.44). The normal form equations are 

dC 2 dU 2 dV 2 

9U 2 ff/ ehh' d( 2 h> s 

-w - fV2 = ~ 9 hTh^ ~ (7 - 8 - 61) 

W + fU2 = ~ 9 hTv^ ~ h T * ' (7 - 8 - 62) 

Note that the effective gravity 

9h = 9^ (7.8.63) 
is much smaller than g. Moreover, the effective wind stresses are 



hf_ s h' 

hence are oppposite to the surface wind stresses. 



* -ji (7.8.64) 



7.8.4 Free modes without wind 

Let us consider free waves in the absence of wind forcing ((r x ) k , (j y ) k ) = 0). It is then 
possible to consider one mode and assume that the other modes to be zero. In particular we 
shall examine (eigen, or natural) modes of sinusoidal wave form 

(Cfc, U k , V k ) oc e ^x+ikyy-iu,t (7.8.65) 
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Barotropic (surface) mode: 

Let us focus on mode 1 alone and assume mode 2 to be absent, i.e., 

U 2 = V 2 = ( 2 = (7.8.66) 

then 

U 2 = -(h'/h)U + U' = 

V 2 = -(h'/h)V + V = (7.8.67) 

C 2 = -(ti/h)( + (i + ti/h)(' = o, 

hence 

^=4 "=4 c '^ c - (7 - 8 - 68) 

Again the depth-averaged velocities in both layers are in phase and are numerically the same. 
The interface displacement is in phase with and smaller than the free surface displacement 
in proportion to the vertical distance above the seabed. 

Moreover, (7.8.54) to (7.8.56) become for the upper layer: 

l(^)-'(^)-^*>§ 



or 



dV dC 
d( h + h' IdU dV\ „ 



It is easy to get the equations for the lower layer. 



dU' dC 

™- f V> = - gh > S A (7.8.75) 
dV dC 

d(' h + ti (dU' dV'\ ,„ m . 
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Consider sinusidodal waves and the equations for the upper layer, 

—iuU — fV = —ghik x ( 



+ 



—iojV — fU = —ghiky( 
h + h',., 



h 



(ik x U + ik y V) = 



For nontrivbil solution th eeigenvalue condition is easily found to be 

UJ 2 = f + g(h + ti)k 2 

which can be written as 



h + h' 
uj\ = ± 



lh + h' 



r\2 
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k\h + h') 



1/2 



(7.8.78) 



(7.8.79) 



This is the dispersion relation for the barotropic (surface wave) mode, and is plotted in 
Figure 7.8.2 



Baroclinic (internal wave) mode 

Now consider the unforced Mode 2 and assume Mode 1 to be absent, 

Ui = Vi = Ci = o 

Since 



ai = 1 



we have, from (7.8.31) to (7.8.33) 



U + U' = 0, V + V = 0. 
eh' 



C 



-C 



(7.8.80) 
(7.8.81) 

(7.8.82) 
(7.8.83) 



h + h' 

Note that the horizontal velocities in the two layers are oppositive in phase, so are the free 
surface and the interface. The last result implies the free surface behaves like a rigid lid; the 
interface oscillates with much greater vigor. 
Using these results, we get 



U 2 

(2 



(7.8.84) 
(7.8.85) 
(7.8.86) 
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Figure 7.8.2: Dispersion relation for the barotropic mode in a two-layered sea of contant 
depth. 

Putting all this in (7.8.34) to (7.8.36), we get 
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Therefore, for the lower layer, 



Similarly, 

Also from (7.8.31), 



dU' 
dV 

~~bT 

h'\ dC 



eh'h h + h' df 
h + h' h dx' 

eh'h d(' 



fV = -g 



+ fU' 



h + h' dx 
geh'h dC 



h J dt \ h J \ dx ' dy 



h + h' dy 

h'\ fdU' dV r 
+ 



= 0. 



(7.8.87) 

(7.8.88) 
(7.8.89) 

(7.8.90) 
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hence, 

dC_ dU^ dV^ _ Q (7 8 91) 

dt dx dy 
We can also write the equations for the upper layer, 

dV rTT geh'h df 

w +fu= hvi (78 - 93) 

and 

d(' _dU _dV 
dt dx dy 

Note the sign changes in the momentum and mass conservation equations. 
For sinusoical waves we have 

—iuU — fV = g'hik x C' 

-iujV -fU = g'hik y C 
-iuC - (ik x U + ikyV) = 

with 

' = *hTW (7 ' 8 ' 96) 

The eigenvalue condition acan be easily found to be 

u 2 = f 2 + g'hk 2 (7.8.96) 

which can be written as 

a V^ =± W? + * 2 ' i2 ) 1/2 <7 - 8 - 97) 

This dispersion relation is similar to Figure 7.8.2. Since the reduced gravity, i svery small. 
The baroclinic mode has a much longer natural period. 
Mode 2 is called the barolinic (internal wave) mode. 



